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A two-dimensional nematic phase of magnetic nanorod^ 

Kostyantyn Slyusarenko,^ Doru Constantin,and Patrick Davidson^ 

Laboratoire de Physique des SoUdes, Universite Paris-Sud, CNRS, UMR 8502, 91405 Orsay, France. 

We report a hybrid mesophase consisting of magnetic nanorods confined between the non-ionic surfactant 
bilayers of a lamellar phase. The magnetic field-induced ordering of the nanorods was measured experimentally 
and modeled by a two-dimensional Onsager theory including the third virial coefficient. The nanorods are 
strongly confined in layers, with no orientational coupling from one layer to the next. At high volume 
concentration they exhibit spontaneous in-plane orientational ordering and form a stack of independent two- 
dimensional nematic systems. This isotropic-nematic transition is first-order. 

PACS numbers: 82.70.Dd, 61.30.-v, 61.05.cf, 78.20.Fm 


I. INTRODUCTION 

Ordered phases in reduced dimensions are fascinat¬ 
ing systems, with physical properties that can be qual¬ 
itatively different from those of their three-dimensional 
counterparts. One of the most (intuitively) straightfor¬ 
ward examples would be the two-dimensional nematic 
phase, with quadrupolar orientational order but no posi¬ 
tional order. This deceivingly simple system has been the 
topic of extensive theoreticalii^ and numerical^^— work, 
with a particular emphasis on the first- or second-order 
character of the isotropic-nematic transition. 

There are however surprisingly few experimental re¬ 
alisations of such phases. The first qualitative observa¬ 
tion involved rigid phospholipid tubules at the air-water 
interface^. Similar results were obtained for Langmuir- 
Blodgett films of nanorods^i^, but the particles could 
only be observed after deposition on a substrate. An¬ 
other strategy relies on inserting DNA molecules within 
stacks of neutral^iiS or chargedii lipid bilayers. These 
disordered mixtures were studied by X-ray scattering and 
the nematic phase was identified indirectly, by analyz¬ 
ing the shape of the DNA interaction peak. It is not 
clear whether these are purely two-dimensional phases 
or whether there is some interlayer coupling. We empha¬ 
size that the nematic order parameter was not measured 
in any of the systems above. 

In this paper, we present a system where rigid mag¬ 
netic nanorods are inserted into a soft lamellar matrix of 
nonionic surfactant. The resulting phase has three key 
advantages: 

• It is easily aligned by thermal treatment, so that 
in the small-angle x-ray scattering (SAXS) images 
we can clearly discriminate between the direction 
of the smectic director and that in the plane of the 
layers. We can thus explore the anisotropy of the 
phase. 
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• The nanorods have magnetic properties, hence 
their orientation can be controlled using an applied 
magnetic field. 

• When the phase is aligned in homeotropic anchor¬ 
ing, we can perform optical birefringence measure¬ 
ments (in particular under magnetic field), which 
yield a completely independent estimate of the ori¬ 
entational order of the particles. 

The formulation of the phasei^ and its structural studyi^ 
have already been presented. Here, we are concerned 
with the orientational order of the particles under mag¬ 
netic field, that we describe using an Onsager-type model 
(up to the third virial coefficient), accounting for the elec¬ 
trostatic effects and for an applied magnetic field. 

We conclusively demonstrate a first-order phase tran¬ 
sition between two-dimensional isotropic and nematic 
phases, settling a long-standing debate in the theoreti¬ 
cal and numerical literatures^— . 


II. MODEL AND ANALYSIS 

The nanorods can be seen as parallelepipeds, with di¬ 
mensions L X W X H and volume V = LWH. They 
bear a permanent magnetic moment /i along their long 
axis and exhibit negative susceptibility anisotropy Ay = 
X|| “ X-L < 0. For simplicity, we will describe them 
in the following as cylinders of length L and diameter 
D = 

The particle orientation is quantified by the distri¬ 
bution /(H), such that a fraction /(H)dH of particles 
have their long axes e within the solid angle element dH 
around direction H. 

Our main result concerns the orientational distribu¬ 
tion of the particles under confinement, but we start the 
discussion with their behavior in aqueous solution, a con¬ 
figuration that we will use as reference. 

In isotropic solution at a volume concentration and 
in the presence of a magnetic field (taken along the polar 
axis, B = Bz), /(H) can be described by an extended 
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Onsager mode l^^i^^ : 


/(^) = j exp 


KB cos 0 + JB^P 2 (cos 9) 

-f 


( 1 ) 


In Q the third and fourth terms on the right-hand 
side describe the second and third virial coefficients of 
the Onsager interaction. The coefficient fc 2 (<p, ‘p') has an 
analytical form, 


k2{tf, ip') 



I sin((p 


ip')\ -I- 2[1 -I- I cos((/7 




with K = and J = 3 ^^^^ , where P 2 is the Legendre 
polynomial of the second order, 7 is the angle between 
two rods oriented along directions fl = {9, (p) and fl' = 
{9', ip'), while 


4D^L 

{D + X){L + X)^ 


( 2 ) 


is the spinodal concentration (above which the isotropic 
state is absolutely unstable), X is an effective length 
describing the electrostatic interaction between rodsi^ 
{X = (In A' -b 7 -b In 2 — 0.5 )/k, with A' the amplitude 
of the electrostatic interaction, 7 Ri 0.577 Euler’s con¬ 
stant, and the Debye screening length), and ^ is a 
normalization constant imposing J f{Q)dfl = 1 . 

Equation o is implicit in /(fl) and must be solved 
numerically. The resulting order parameter S is used to 
fit the experimental data. 

When the particles are confined in the lamellar phase 
(with interlayer spacing d), we consider that within the 
same layer they have the same effective interaction as 
in solution so we model them using the 2D analogue of 
Equation o, but this time the polar axis is along the 
director z\\n and the magnetic field is applied in the plane 
of the layers BTn : 


where A 


L + X 
D + X 


is the effective aspect ratio, and 


the coefficient k'i{ip,ip',ip") is calculated numerically^. 
Note that Onsager’s argument for neglecting higher-order 
virial terms does not apply in 2 Did^, especially for rods 
of moderate aspect ratio (such as those used in our ex¬ 
periments). We therefore kept the third term in the virial 
expansion. 

As presented so far, the 2D model considers that parti¬ 
cles in the same plane interact by the effective hard-core 
model measured in solution, while particles in different 
layers do not interact at all. This latter feature is indis¬ 
pensable for having a true 2D phase, so it should be thor¬ 
oughly checked. That is why we added in ([S]) a possible 
orientational couplingi^ between the nanorods in neigh¬ 
boring layers, described by the last term (with M the 
coupling amplitude). Below, we discuss possible physi¬ 
cal origins for such a coupling, estimate its strength and 
show that it is negligible. 


A. X-ray scattering 

The X-ray scattering form factor of an individual par¬ 
ticle is: 


fi^) = ^ exp 
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where the in-plane angle p is measured with respect to 
the orientation of the magnetic field and 7 , 


+ X){D + X) 

^ ~ ttD^L ’ 

is the effective surface fraction. By using D, we implicitly 
consider that the particles can rotate freely along their 
long axis, as in the isotropic solution. If they were to 
adopt a particular configuration, for instance with the 
largest face (L x IT) in contact with the surfactant bi¬ 
layer, the values of X and rj would be slightly changed. 
However, we have no experimental evidence for or against 
this assertion and, for simplicity’s sake, we conserve the 
“cylindrical” model used above. 


iPt a T IP ■ f qL cos I3\ f Ji{qD/2 sin I3)\ 

f (,, ri,L,D) = F, »„c j ( j 

( 4 ) 

where /3 is the angle between the long axis of the rod 
e and the scattering vector q: cos/3 = cos0sin 7 /; -b 
sin 0 cos t/i cos <73 (see also Figure [1]). 

For a monodisperse population, the scattering signal is 
given by 2 £: 

p7Tj2 

I{q,'tlj) = 2 dp /Ms(0)E^(g,/3)sin(0)d0 (5) 
Jo Jo 

while in the general (polydisperse) case Equation (O 
must be averaged over a radius distribution g{D) (the 
polydispersity in the length L does not contribute over 
the available g-range.) In ([5]) we approximate the orien¬ 
tational distribution of the rods by the analytical Maier- 
Saupe form: 

/„,(n) = ( 8 ) 

where Z is a normalization constant and fitted the data 
with only one parameter m. The order parameter S can 
be expressed analytically as a function of m. 
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FIG. 1. Experimental geometry. The scattering object is a 
rod with orientation e given by the polar angle 6 (measured 
with respect to the nematic director n|| 2 ) and the azimuthal 
angle ip. The incident beam is along the y axis. The detector 
records in the {x, z) plane (normal to the incident beam) the 
scattered intensity I{q, ip) given by Equation ([5| 

B. Optical birefringence 

The birefringence of the samples can generally be ex¬ 
pressed as: 

An = Ansi,tS{B)4> (7) 

where the field-dependent order parameter S{B) is the 
uniaxial one S for the case of water suspensions and the 
in-plane component |P of the biaxial order when the 
particles are confined within the lamellar matrix. Arisat 
is the specific birefringence, depending on the geometry 
of the particles and on their dielectric permittivity tensor. 


III. RESULTS AND DISCUSSION 

From electron microscopy and X-ray diffraction mea¬ 
surements, the average particle dimensions are: 

L = 315 nm x W = 38 nm x H = 18 nm, ( 8 ) 

and the relative polydispersity is of the order of 0.3^^. 
They correspond to cylinders with a diameter D = 
= 29.5 nm. 



FIG. 2. Order parameter S{B) (determined from the SAXS 
images) for the aqueous suspensions of nanorods at various 
concentrations 0 (symbols). The curves are fits derived from 
Equation m with parameter values 


The order parameter is then computed over the distri¬ 
bution f{9) determined from the complete model ([T]) as 
a function of B and cp, yielding estimates for the material 
constants K and J. 

The spinodal concentration (p* = 7.1% was estimated 
from the independent measurement of the phase diagram 
of the aqueous suspension in zero field: (p* ~ 0.950^, 
where is the volume concentration of the nematic 
phase at coexistence (in aqueous solution). From ([2]) we 
then obtain X = 72 nm. 

We also measured the birefringence An{B) of the so¬ 
lutions (Figure |3]) and fitted it to Equation ([7]), which 
involves the additional constant Arisat- 

Finally, we obtain the material constants: 


A: = 8.2±0.5T-i (9a) 

J = -13.8±IT-2 (9b) 

Arisat = 0.8 (9c) 

P* = 7.1% ^ X = 72 nm (9d) 


A. In solution 

Under a magnetic field B, aqueous solutions of 
nanoparticles become uniaxial, with a director n||B and 
an order parameter S that can be positive or negative, 
depending on the field amplitude^^. 

From the SAXS images we determined the order pa¬ 
rameter S{B) as a function of the field (Figure [2]) by 
fitting the data to Equation ([5]) and using for the orien¬ 
tation distribution the Maier-Saupe form ([ 6 |). The poly¬ 
dispersity 5 (Z?) is estimated from the scattering spectrum 
of a dilute suspension in the absence of the field. 


that we use alongside the particle dimensions ([ 8 ]) to study 
the hybrid system (particles inserted in the lamellar ma¬ 
trix). 


B. In the lamellar phase 

The SAXS measurements can be performed in either 
homeotropic or planar anchoring, giving access to both 
S and P, but the optical birefringence can only be de¬ 
termined in the homeotropic configuration (sensitive to 

P)- 
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FIG. 3. Birefringence An{B) of the aqueous suspensions of 
nanorods at various concentrations (p (symbols). The curves 
are fits derived from Equations © and 0 with parameter 
values 0 


In the planar configuration, with B = 0, for (p = 0.5%, 
1.5%, and 3.5% we obtained S = —0.44, —0.46, and 
—0.47, respectively, very close to the perfect confinement 
limit S = —0.5. We conclude that, even in the absence 
of the field, the particles are practically contained within 
the plane of the layers (and the optical axis is along n). 
If we now apply a field B±n the system becomes biaxial, 
with a second order parameter P. 

In the following we treat the system as purely two- 
dimensional, with S = —0.5. The particle population is 
then completely described by the orientation distribution 
/((^), determined from the model (0 as a function of B, 
M, and (p. The material constants are those measured 
in aqueous solution (0, the periodicity of the lamellar 
phase d = 45 nm,™ and M is the sole fitting parameter. 


1. Isotropic two-dimensional phase 

We describe the birefringence An{B) using Equa¬ 
tion 0 (with S{B) replaced by (2/3)P(B)) coupled 
with 0. The simultaneous fit of the three curves yields 
M = {—3±3)kBT see Figure 0 and the resulting or¬ 
der parameter is compared to that measured by SAXS in 
Figure [5] The M values for the individual fits at 1.5 and 
3.5 % are shown in Figure© At 0.5 % the uncertainty 
on M is very large, so we excluded this point from the 
analysis. 


2. Nematic two-dimensional phase 

The most concentrated system, with (p = 8 vol % in the 
lamellar phase is a biaxial nematic, with a spontaneous 



B (ml) 


FIG. 4. Birefringence An{B) of the suspensions of nanorods 
in lamellar phase at various concentrations p (symbols). The 
curves are fits derived from Equations 0 (with S{B) replaced 
by {2/3)P{B)) and 0 with the material parameters in 0 , 
d = 45 nm, and M = —3 ksT. 



FIG. 5. Second order parameter P{B) (determined from the 
SAXS images) of the suspensions of nanorods in lamellar 
phase at various concentrations p (symbols). The curves are 
fits derived from Equation @ with the material parameters 
in 0, d = 45 nm, and M = —SksT. 

second order parameter P. It can thus be seen as a stack 
of two-dimensional nematic layers. To measure P we ap¬ 
plied a small magnetic field of 30 mT in the plane of the 
layers to orient all nematic domains along the field, and 
measured the dependence of the x-ray scattering intensity 
/ on pj. Fitting I(ip) by the Leadbetter method^i, which 
takes into account the interaction between nanorods (and 
is more appropriate at high concentration than the model 
we obtained P = 0.51. 

Both in aqueous solution and in the lamellar phase 




















5 


the isotropic-nematic transition is first-order, as the two 
phases can coexist over a certain concentration range. 
By preparing several samples within this domain and es¬ 
timating the fraction occupied by each phase we deter¬ 
mine the nematic concentration at coexistence (j>^ (in 
the aqueous solution) and (pN (for the two-dimensional 
nematic). In our system, these two distinct parameters 
happen to have the same numerical value, 7.5±0.5 vol %. 
In the lamellar phase, the isotropic-nematic coexistence 
range is approximately p = 4.5 — 7.5 vol %, correspond¬ 
ing to r] = 0.38 — 0.63i^. A sample at coexistence, with 
4> = 6.6 vol %, i.e. r] — 0.54, is shown in FigurelHl The in¬ 
terface is not sharply defined, mainly due to the presence 
of the lamellar defects which create arbitrarily shaped 
domains of the two phases. 



FIG. 6 . Sample at coexistence {p = 6.6 vol %, 77 = 0.54). The 
dashed line (with a greyed uncertain area) delimits the 2D 
isotropic phase I/La (upper left) from the 2D nematic phase 
N/La (lower right). The width of the capillary is 1 mm. 

The two-dimensional model ([3]) does indeed predict a 
first-order isotropic-nematic transition, with a (pN that 
corresponds to the experimental one ((/>™^ = 7.5 ± 
0.5 vol %) for M = (0.5 ± 0.5) LbT. For our most con¬ 
centrated sample, with p = 8vol%, the same model 
reproduces the order parameter P = 0.51 for M = 
(0.07 ± 0.01) ksT. These two values correspond to the 
high-concentration points in Figure [T] 


3. Coupling 

The coupling coefficient M is not significantly different 
from zero: at low concentration (in the two-dimensional 
isotropic phase) M = (—3±3) ksT from the birefringence 
data and at high concentration (in the two-dimensional 
nematic phase) M = (0.5 ± 0.5) ksT from the tran¬ 
sition concentration and M = (0.07 ± 0.01) ksT from 
P{p = 0.08). We conclude that, within the experimental 
precision, the orientation of particles in different layers is 
uncoupled, see Figured Is this result in agreement with 
the theoretical estimations? 



FIG. 7. Best estimate of M for various concentrations p. For 
p = 1.5 and 3.5 % we used the birefringence data (Figure 
0 . The two high-concentration points are obtained from the 
transition concentration and from P[p = 0.08) (see text). 


The main difference between the particles in lamel¬ 
lar phase and those in isotropic solution is that the for¬ 
mer are constrained to lie in parallel planes, and thus 
can more easily become parallel. Magnetic and electro¬ 
static interaction might then be stronger than in three 
dimensions. We estimate their amplitude in the two- 
dimensional nematic phase, at </> = 8 vol %: 

The magnetic interaction between the permanent 
dipole moments of the particles IFdd favors an antiparal¬ 
lel orientation and its magnitude is: 

ITdd _ _ ^loK'^ksTp _ ^ ^ g 

ksT Airr^ksT AttLDH ’ 

where r is an average interparticle distance. 

The electrostatic twist energy2^ per molecule Wtw 
(which favors a perpendicular orientation) is given by2^: 

y J (10) 

where /itw = The coefficient ktwi^P,^^') is: 


kt^{(p,ip') = 

— a\ sin((^ — 





[-|sin((^ 


ip') \ In I sin(i75 




with 

/ / I sin(tp - ‘P')\ In I sin((/? - ip')\d(pdip' _ 

“ / / I sin((^ - ip')\dipd^' ^ 

The main uncertainty is in evaluating the Debye length: 
K,~^ = 2 — 14 nm, where the lower value was measured in a 
different goethite batchii and the higher one corresponds 
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to an ionic strength of 1 mM (the minimum possible at 
pH = 3). We then have Wtw ~ (4 — 25) • 10 “^A:bT. Both 
effects are therefore negligible compared to the thermal 
energy, in agreement with our experimental findings. 


IV. CONCLUSION 

We formulated stable anisotropic materials where 
goethite nanorods are confined between the bilayers of 
a soft lamellar phase. The nanorods form a nematic sys- 


temlwith director along the director of the lamellar phase 
and a negative order parameter — U.4t) (to be com- 
pared with S = —0.5 for perfect confinement). 

At low volume concentration (j) < 7.5 vol %, the system 
is i niiaxiaf in the absence of a magnetic field and forms a 
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ORDER PARAMETER 
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the distribution /(fl)) and can be written in the princi¬ 
pal axis frame as (Palffy-Muhoray and Hoatson—, Eq. 
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(-\{S-P) 0 0\ 

Qa/3= 0 -i(5 + P) 0 (11) 
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In spherical coordinates, with the polar axis along z, 
the values of P and S are related to the distribution 
function f{0, ip) via: 



f{9, p) sin 9 cos^ 9 d9Ap — ^ 
/(0, p) sin^ 9 cos(2(^) d9d.p 


( 12 ) 


For a uniaxial system the distribution / depends only on 
the polar angle, /(H) = f{9), P = 0 and S is defined 
via the simplified formula: S = ^ f f{9) sin^cos^ 9d9 — 
y In the “complete confinement” case (relevant for the 
lamellar system), /(H) = 6{9 — n/2)f{p), S = and 
P reduces to: P = § f f{p) cos(2i^) dp. 







